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Abstract 

We consider the late-time evolution of gravitational and electromagnetic per- 
turbations in realistic rotating Kerr spacetimes. We give a detailed analysis 



^r^ ' of the mode-coupling phenomena in rotating gravitational collapse. A con- 

^^ ' sequence of this phenomena is that the late-time tail is dominated by modes 

o 



which, in general, may have an angular distribution different from the original 
one. In addition, we show that different types of fields have different decaying 
rates. This result turns over the traditional belief (which has been widely 
accepted during the last three decades) that the late-time tail of gravitational 



/\ • collapse is universal. 

I. INTRODUCTION 

The no-hair conjecture, introduced by Wheeler in the early 1970s |I[], states that the 
external field of a black hole relaxes to a Kerr-Newman field characterized solely by the 
black-hole mass, charge and angular momentum. 

Price was the first to analyze the mechanism by which the spacetime outside a 
(nearly spherical) star divests itself of all radiative multipole moments, and leaves behind a 
Schwarzschild black hole; it was demonstrated that all radiative perturbations decay asymp- 
totically as an inverse power of time, the power indices equal 2/ + 3 (in absolute value), 



where / is the inultipole order of the perturbation. These inverse power-law tails are a direct 
physical consequence of the backscattering of waves off the effective curvature potential at 
asymptotically far regions 0,0. Leaver |^ demonstrated that the late-time tail can be as- 
sociated mathematically with the existence of a branch cut in the Green's function for the 
wave propagation problem. 

The analysis of Price has been extended by many authors. We shall not attempt to 
review the numerous works that have been written addressing the problem of the late-time 
evolution of gravitational collapse. For a partial list of references, see e.g., |5HT6|. 



The above mentioned analyses were restricted, however, to spherically symmetric back- 
grounds. It is well-known, however, that realistic stellar objects generally rotate about their 
axis, and are therefore not spherical. Thus, the nature of the physical process of stellar 
core collapse to form a black hole is essentially non-spheric, and an astrophysically realistic 
model must take into account the angular momentum of the background geometry. 

The corresponding problem of wave dynamics in realistic rotating Kerr spacetimes is 
much more complicated due to the lack of spherical symmetry. A first progress has been 
achieved only recently [|T7|-|2^. Evidently, the most interesting situation from a physical 
point of view is the dynamics of gravitational waves in rotating Kerr spacetimes. Recently, 
we have begun an analytic study of this fascinating problem ||2^ . This was done by analyzing 



the asymptotic late-time solutions of Teukolsky's master equation pj,|25|, which governs the 
evolution of massless perturbations fields in Kerr spacetimes. In this paper we give a detailed 
analysis of the problem. In particular, we give a full account to the phenomena of mode 
coupling in rotating spacetimes (this phenomena has been observed in numerical solutions 



of Teukolsky's equation |jT^,|T8[) 

The plan of the paper is as follows. In Sec. |I| we give a short description of the 
physical system and summarize the main analytical results presented in Ref. [^. In Sec. 
pn] we discuss the effects of rotation and the mathematical tools needed for the physical 
analysis are derived. In Sec. ^ we analyze the active coupling of different gravitational and 
electromagnetic modes during a rotating gravitational collapse, with pure initial data. In 



Sec. we consider the late-time evolution of realistic rotating gravitational collapse, with 
generic initial data. We conclude in Sec. [V| with a summary of our analytical results and 
their physical implications. 

II. REVIEW OF RECENT ANALYTICAL RESULTS 



The dynamics of massless perturbations outside a realistic rotating Kerr black hole is 
governed by Teukolsky's master equation [^,^: 
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where M and a are the mass and angular-momentum per unit-mass of the black hole, and 
A = r^ — 2Mr + a^. (We use gravitational units in which G = c = 1). The parameter 
s is called the spin-weight of the field. For gravitational perturbations s = ±2, while for 
electromagnetic perturbations s = ±1. The field quantities ip which satisfy Teukolsky's 
equation are given in [^]. 

Resolving the field in the form 



^ = A-^/2(r2 + a2)-i/2 Y. ^"e 
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where m is the azimuthal number, one obtains a wave-equation for each value of m (we 
suppress the index m) 
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where the tortoise radial coordinate y is defined by dy = ^ ^" dr. The coefficients Bi{r, 
are given by 
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[The explicit expression of B^{r,6) is not important for the analysis]. 
The time-evolution of a wave-field described by Eq. (^) is given by 

m{z, t) = 2n [ r \b^{z') \G{z, z'- t)^t{z', 0) + Gt{z, z'- t)m{z\ 0) 



+ 



B2{z')G{z, z'- t)^{z\ 0) \?,me'de'dy' 
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for t > 0, where z stands for (y, 9). The (retarded) Green's function G{z, z'\ t) is defined by 
DG{z, z'; t) = 5{t)S{y - y')5{9 - 9')/27ism9, with G = for t < 0. We express the Green's 
function in terms of the the Fourier transform Gi{y, y'; u) 
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where c is some positive constant and /q = max(|m|, |s|). The functions sSp{9,auj) are the 
spin-weighted spheroidal harmonics which are solutions to the angular equation 
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For the au = case, the eigenf unctions sSl"{9,au) reduce to the spin-weighted spherical 
harmonics sy/'"(6', 0) = sS]^{9)e'^^'^, and the separation constants sA^i^au) are simply gA]^ = 
il-s){l + s + l) H. 

The Fourier transform is analytic in the upper half if-plane and it satisfies the equation 
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wherein = {r^+a?)ixi—am, X = A+a^uj^—2amuj, and H = s{r—M)/{r'^+a )+rA/{r^+a^) 



Define two auxiliary functions \l/i and \l/2 as solutions to the homogeneous equation 
D{uj)'^i=D{uj)'$2 = 0, with the physical boundary conditions of purely ingoing waves cross- 
ing the event horizon, and purely outgoing waves at spatial infinity, respectively. In terms of 
\^i and ^2, and henceforth assuming y' < y, Gi{y, y'; u) = —'^i{y', uj)'^2{y, ^)/W{uj), where 
we have used the Wronskian relation W{uj) = W^ (^i, ^2) = ^i^2,y - ^2^i,y 

It is well-known that the late-time behaviour of massless perturbations fields is de- 
termined by the backscattering from asymptotically far regions P,^. Thus, the late- 
time behaviour is dominated by the /oiy-frequencies contribution to the Green's function, 
for only low frequencies will be backscattered by the small effective curvature potential 
(at r ^ M). Therefore, a small-uj approximation [or equivalently, a large-r approxi- 
mation of Eq. (^)] is sufficient in order to study the asymptotic late-time behaviour of 
the fields [|12|. With this approximation, the two basic solutions required in order to 
build the Fourier transform are \&i = r'"''^e*'^''M(/ + s + 1 — 2iujM,2l + 2,—2iujr), and 
^2 = r^+^e''^''U{l + S + 1- 2iwM, 21 + 2, -2iijr), where M{a, b, z) and U{a, b, z) are the two 
standard solutions to the confluent hypergeometric equation [^. Then 

W{^i, ^2) = t{-iy+\2l + l)!(2cu)-(2'+i)/(/ + s)\ . (10) 

In order to calculate G{z,z';t) using Eq. (|^), one may close the contour of integration 
into the lower half of the complex frequency plane. Then, one identifies three distinct con- 
tributions to G{z, z'; t) 0] : Prompt contribution, quasinormal modes, and tail contribution. 
The late-time tail is associated with the existence of a branch cut (in \I'2) in the complex fre- 
quency plane [Q (usually placed along the negative imaginary u— axis). A little arithmetic 
leads to p3| 

"^2(l/,cue2-) ^2iy,uj) 
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III. ROTATION EFFECTS THE COUPLING OF DIFFERENT MODES 

The rotational dragging of reference frames, caused by the rotation of the black hole (or 
star) produces an active coupling between modes of different I (but the same m). Mathe- 
matically, it is the 6'-dependence of the spin-weighted spheroidal wave functions sSl^{6,auj) 
and of the coefficients Bi{r, 6) and B2{r, 6) which is responsible for the interaction between 
different modes; no coupling occurs in the non-rotating (a = 0) case. 

The angular equation (§) is amenable to a perturbation treatment for small aw p8| , p9 



we write it in the form {L^ + L^)sSl" = —gAf^sSp, where L^{0) is the uj- independent pait of 
Eq. (D, and 

L^{6,auj) = {aujYcos'^9 — 2auscos9 , (13) 

and we use the spin-weighted spherical functions sY{^ as a representation. They satisfy 
L'^sYi = —sA] sYi with sA} = (/ — s)(/ -|- s -|- 1) (we suppress the index m on gAi and sYi)- 



For small acu a standard perturbation theory yields (see, for example, |^ 



Si{e, au) = g Cikiauy'-'^KYkie) , (14) 



sJl 



k=lo 



where, to leading order in aoo, the coefficients Cik{auj) are uj— independent pT| , |29| . Equation 
(p!^) implies that the black-hole rotation mixes (and ignites) different spin-weighted spherical 
harmonics. 

The coefficients -Bi(r, 6) and B2{r, 9) appearing in the time-evolution equation (|^) depend 
explicitly on the angular variable 9 through the rotation of the black hole (no such depen- 
dence exist in the non-rotating a = case). Therefore, in order to elucidate the coupling 
between different modes we should evaluate the integrals {slm\skm) , {slm\sm'^9\skm) , and 
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(s/m|cos^|sA;m), where {slm\F{e)\skm) = J sY{^*F{e)sY,^dn [see Eqs. (§ and (D for the 
definition of the Bi{r,6) coefficients]. In addition, the values of the coefficients Cik depend 
on the integrals pT| , p9| {sl'm\cos'^6\skm) and {slmlcosOlskm) [see Eq. (|I3|) for the definition 
of the perturbation term L^{6, auj), which is responsible for the mixing of modes in rotating 
backgrounds] . 

The spin-weighted spherical harmonics are related to the rotation matrix elements of 
quantum mechanics |3l|. Hence, standard formulae are available for integrating the prod- 
uct of three such functions (these are given in terms of the Clebsch-Gordan coefficients 
28| , |2l|j29| ) ■ In particular, the integrals (s/0|sin^^|sA;0) and (s/0|cos^^|sA;0) vanish unless 



/ = k,k ±2, while the integral (s/0|cos^|s/cO) vanishes unless I = k ± 1. For non-axially 
symmetric {m 7^ 0) modes, {slm\sm'^9\sk'm) j^ for I = k, k ± 1, k ± 2 (the same holds 
for the integral (s/0|cos^6'|s/i;0)), and (s/?Ti|cos^|s/i;?Ti) 7^ for I = k,k ±1 (all other matrix 
elements vanish). Note also that the complex coefficient B2 couples the real and imaginary 
parts of \E'"^. 

We are now in a position to evaluate the late-time evolution of realistic rotating grav- 
itational collapse. We shall consider two kinds of initial data: Pure initial data, which 
corresponds to the assumption that the initial angular distribution is characterized by a 
pure spin-weighted spherical harmonic function sV/™, and generic initial where the initial 
pulse consists of all allowed modes (all spherical harmonics functions with / > /q) • 

IV. PURE INITIAL DATA 

A. Asymptotic behaviour at timelike infinity 

As explained, the late-time behaviour of the fields should follow from the /oty-frequency 
contribution to the Green's function. Actually, it is easy to verify that the effective con- 
tribution to the integral in Eq. (0) should come from \uj\=0{l/t). Thus, we may use 
the Iti^lr <^ 1 limit of \E'i(r, cj) in order to obtain the asymptotic behaviour of the fields at 



timelike infinity (where y,y' ^ t). Using Eq. 13.5.5 of [|2^ one finds \l/i(r, cu) ~ Ar^^^. 



Substituting this in Eq. ([12|), and using the representation Eq. ([1^) for the spin- weighted 
spheroidal wave functions gSi, togather with the cited properties of the angular integrals (of 
the form {slm\F{6)\skm)), we find that the asymptotic late-time behaviour of the / mode 
(where / > /q) is dominated by the following effective Green's function: 

, ._ f M(-l)(^'+^+^-''-^-)/^2^^+^(fc + .)!(A: -.)!(/* + / + 2 - q)\ 

' ^ ' ' ^ " h Ai'^k + iw ^^^ ^ 

where q = min(/* — Zq, 2). Here, L = /* — g for / > /* — 1 modes, and L = / for / < /* — 2 
modes. Thus, the late-time behaviour of the gravitational and electromagnetic fields at the 
asymptotic region of timelike infinity i^ is dominated by the lowest allowed mode, i.e., by 
the / = Iq mode. The corresponding damping exponent is — (/* + /q + 3 — g). 

B. Asymptotic behaviour at future null infinity 

We further consider the behaviour of the fields at the asymptotic region of future null 
infinity scri^. It is easy to verify that for this case the effective frequencies contribut- 
ing to the integral in Eq. (O) are of order Oilju). Thus, ior y — y' <^ t <^ 2y — y' 



one may use the \uj\y' <^ 1 asymptotic limit of ^i{y',u!) and the M <^ \uj\ ^ '^ y 
[Imuj < 0) asymptotic limit of ^i{y,u!). Thus, \l'i(|/', cj) ~ Ay''-^^, and ^i{y,uj) ~ 
e'^y{2l + l)!e-^'^('+*+^)/2(^2w)-('+^+^)2/-V(/ - s)\, where we have used Eqs. 13.5.5 and 13.5.1 
of ||2^, respectively. Substituting this in Eq. (p!2D , and using the representation Eq. (p!4D 



for the spin- weighted spheroidal wave functions gSi, togather with the cited properties of 
the angular integrals, one finds that the behaviour of the / mode (where I > Iq) along 
the asymptotic region of null infinity scri^ is dominated by the following effective Green's 
functions: 



G?iz,z';t)= E -^^^-^ ^oi \: y 



V 



k=i*-,. 7r(2A: + l)! 

CMsYi{e),Y:{ey''u-^'-^+'^ , (16) 
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for I > I* — 1 modes, where gi = min(/* — Iq, 2) and q2 = niin(/ — /*, 2), and 

_ M(-l)('-+'-^-)/^2^(/ + .)!(/- -.-!)! 
^i[z,z,t)- vr(2/ + l)! y "" 

Cu^^2sYi{9)sY;^,{e')a'*-'-\-^''-^^ , (17) 

for Z < /* — 2 modes. The dominant modes at null infinity and the corresponding damping 
exponents are given in Table |. 

C. Asymptotic behaviour along the black- hole outer horizon 

The asymptotic solution to the homogeneous equation D{u!)'^i{y,u!) = at the black- 
hole outer horizon H+ {y -^ -oo) is [|5[ ^i{y,uj) = C{uj)A-'/'^e''^'^'"''^+^y, where u;+ = 



a/(2Afr+) [r^ = M + (M^ — a^)^/^ is the location of the black-hole outer horizon]. In 
addition, we use ^i{y',u!) ~ Ay''-^^. Regularity of the solution requires C to be an analytic 
function of u. We thus expand C{uj) = Cq + Ciuo + ■ ■ ■ for small uo (as already explained, 
the late-time behaviour of the field is dominated by the /ow>frequency contribution to the 
Green's function). 

Substituting this into Eq. ([l2l), and using the representation Eq. (|lj) for the spin- 
weighted spheroidal wave functions sSu we find that the asymptotic behaviour of the / mode 
(where I > Iq) along the black-hole outer horizon H^ is dominated by the following effective 
Green's function: 

^ f M(-l)(^-+^+^-^-^-)/^2^^+HA: + s)\{k~s)\{l* + l + 2-q)\ ^ _,„^ „^, 
^' ^^' ^ ' ^^ - ,t;„ ' ' 7r[(2fc + l)!p ^ y 



where g,p and L are defined as before, s^k are constants, and 6 = generically, except for 



the unique case m = with s > 0, in which 6=1 ||3^. Hence, the late-time behaviour of 
the gravitational and electromagnetic fields along the black-hole outer horizon is dominated 
by the lowest allowed mode, i.e., by the / = Iq mode. The corresponding damping exponent 

is -{l* + lo + 3-q + b). 



V. GENERIC INITIAL DATA 

So far we have assumed that the initial pulse is made of pure data, characterized by one 
particular spherical harmonic function ^l^r. In this section we consider the generic case. 
That is, we assume that the initial pulse consists of all the allowed (/ > /q) modes (see also 
the most recent analysis of Barack [Q). 

The analysis here is very similar to the one presented in Sec. JV]: Using Eq. (jT^, 
togather with the appropriate asymptotic forms of '^i{y,uj) and '^i{y',uj) (as given in Sec. 
[rV| for the various asymptotic regions), and the representation Eq. (|I^) for the spheroidal 
wave functions, we find that the asymptotic late-time behaviour of the / mode (where / > Iq) 
is dominated by the following effective Green's functions: 

Gf{z, z'- 1) = MF,{yy't+\YmsY,l{d')a'-'H~^'+'-^'^ , (19) 

at timelike infinity i+, where Fi = Fi{l,lo,m,s) = (_i)a+'o+2^+2)/222io+i(/ + Iq + 2)!(/o + 
s)!(/o-s)!G„,/7^[(2/o + l)!]^ 

Gf{z, z'; t)=Yl MF2y"+'v-^Yiie)sY:{ey-'u-^^-^+'^ , (20) 

k=lo 

at future null infinity scri+, where F2 = F2{l,k,m,s) = (_i)0+'«+2s+2)/22A:(^/. _I_ ^yj^i — ^ + 
l)\Cki/iT{2k + l)\,and 

Gf{z, z'; t) = ,r;MFi A-^/2^"«+^y5(^),Y;*(^>'-'"e^'"'^+^i;-('+'°+3+^) , (21) 

at the black-hole outer horizon if+, where s^'i are constants 

VI. SUMMARY AND PHYSICAL IMPLICATIONS 

We have analyzed the dynamics of gravitational (physically, the most interesting case) 
and electromagnetic fields in realistic rotating black-hole spacetimes. The main results and 
their physical implications are as follows: 
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(1) We have shown that the late-time evolution of realistic rotating gravitational collapse 
is characterized by inverse power-law decaying tails at the three asymptotic regions: timelike 
infinity z+, future null infinity scri+, and the black-hole outer-horizon if+ (where the power- 
law behaviour is multiplied by an oscillatory term, caused by the dragging of reference 
frames at the event horizon). The relaxation of the fields is in accordance with the no-hair 
conjecture [|l|. This work reveals the dynamical physical mechanism behind this conjecture 
in the context of rotating gravitational collapse. 

The dominant modes at asymptotic late-times and the values of the corresponding damp- 
ing exponents are summarized in Table | (for pure initial data) and Table (for generic 
initial data). For reference we also include in Table |T| the results for the scalar field toy 



model with pure initial data (the s = case) 0,0] (the results for generic initial data 
coincide with those of gravitational and electromagnetic perturbations). In these tables, / is 
the multipole order of the perturbation, Zq = max(|r7i|, |s|), and /* is the initial mode of the 
perturbation (for pure initial data). For the scalar field case (s = 0), we have p = if / — |m| 
is even, and p = 1 otherwise. Note that for pure initial data, the pulse with /* = /q, ^o + 1 
differs from initial data with Iq + 2 < I* . This is caused by the fact that the Iq mode is not 
ignited (not coupled) to modes with smaller values of /. 

The somewhat different character of the scalar field case can be traced back to Eq. (^ 



for B2{r,6) and Eq. ([131) for L^{6,auj)] it turns out that B2 is 6-independent in the s = 
case, and thus this term cannot couple different modes. To this we should add the fact that 
for the scalar field case, L^{6,auj) is proportional to (ac<j)^ (the term proportional to aws 
vanishes), and thus the coefficients Cik in Eq. (|^) vanish if \l — k\ is odd pT| . 

The damping exponents for generic initial data derived in this paper agree with those 
derived most recently by Barack [Q using an independent analysis. Note, however, that 
Barack's analysis cannot yield the values of the damping exponents for pure initial data. 

(2) The unique and important feature of rotating gravitational collapse (besides the 
oscillatory behaviour along the black-hole horizon) is the active coupling of different modes. 
Physically, this phenomena is caused by the dragging of reference frames, due to the black- 
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hole (or star's) rotation (this phenomena is absent in the non-rotating a = case). As 
a consequence, the late-time evolution of realistic rotating gravitational collapse has an 
angular distribution which is generically different from the original angular distribution (in 
the initial pulse). 

(3) We emphasize that the power indices at a fixed radius in rotating Kerr spacetimes 
(/ + /o + 3 for generic initial data and I* + I + 3 — q for pure initial data) are generically 
sma//er than the corresponding power indices (the well-known 2/+3) in spherically symmetnc 
Schwarzschild spacetimes. (For generic initial data there is an equality only for the lowest 
allowed mode I = Iq, while for pure initial data there is an equality only for the / = /q mode 
provided it characterizes the initial pulse). This implies a slower decay of perturbations in 
rotating Kerr spacetimes. Stated in a more pictorial way, a rotating Kerr black hole becomes 
'bald' slower than a spherically-symmetric Schwarzschild black hole. From Eq. (|1^) it is easy 
to see that the time scale tc at which the late-time tail of rotating gravitational collapse is 
considerably different from the corresponding tail of non rotating collapse (for / > Iq modes) 
is tc = yy'/a, where y' is roughly the average location of the initial pulse. 

(4) It has been widely accepted that the late-time tail of gravitational collapse is uni- 
versal in the sense that it is independent of the type of the massless field considered (e.g., 
scalar, neutrino, electromagnetic, and gravitational). This belief was based on spherically 
symmetric analyses. Our analysis, however, turns over this point of view. In particular, the 
power indices / + /q + 3 at a fixed radius found in our analysis are generically different from 



those obtained in the scalar field toy-model pT| , p2| / + \m\ + p + 3 (where p = ii I — \m 



is even, and p = 1 otherwise). Thus, different types of fields have different decaying-rates. 
This is a rather surprising conclusion, which has been overlooked in the last three decades ! 
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TABLES 
TABLE L Dominant modes and asymptotic damping exponents for gravitational and electro- 
magnetic fields - pure initial data. Iq = max(|m|, \s\), and I* is the initial mode of the perturbation. 



asymptotic region 


I* 


dominant mode(s) 


damping exponent 


timelike infinity 


lo<l* <lo + i 


^0 


-(2/0 + 3) 




h + 2<r 


^0 


-(r + /o + i) 


null infinity 


k<l* <k + l 


^0 


-{lo-s + 2) 




lo + 2<l* 


lo<l<l*-2 


-il*-s) 


outer horizon 


lo<l* <lo + i 


lo 


-(2/0 + 3 + 6) 




lo + 2<l* 


k 


-{l* + lo + l + b) 



TABLE IL Dominant modes and asymptotic damping exponents - generic initial data, 
asymptotic region dominant mode damping exponent 

timelike infinity Zo — (2/o + 3) 

null infinity Zq —{Iq — s + 2) 

outer horizon Zq — (2/o + 3 + 6) 



TABLE in. Dominant modes and asymptotic damping exponents for scalar fields - pure initial 
data, p = if / — |m| is even, and p = 1 otherwise. 



asymptotic region 


I* 


dominant mode(s) 


damping exponent 


timelike infinity 


lo<l* <lo + l 


r 


-(2/* + 3) 




h + 2<r 


lo+P 


-{l* + lo+p + l) 


null infinity 


Iq<1* <lo + l 


I* 


-(r + 2) 




h + 2<r 


lo+P<l<l* -2 


-/* 


outer horizon 


lo<l* <lo + l 


I* 


-(2/* + 3) 




lo + 2<l* 


lo+P 


-{l* + lo+p+l) 
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